AD-A172  779  AN  EXTENSION  OF  THE  k'REISS  MATRIX  THEOREMCU)  UISCONSIN  1/1 
UUN I V-MAD I  SON  MATHEMATICS  RESEARCH  CENTER 
B  A  UADE  ET  AL  AUG  8b  MRC-TSR-2951  DAA229-80-C-0041 

F/G  12/1 


UNCLASSIFIED 


NL 


1.0  Bf* 


M  12.5 


IM  lii  Hi  1.6 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS  1963  A 


MRC  Technical  Summary  Report  #2951 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


August  1986 


(Received  August  13,  1986) 


DTIC 

SELECTE| 
OCT  8  1986 

B 


Approved  for  public  release 
Distribution  unlimited 


onsored  by 


Army  Research  Office 
Box  12211 


Research  Triangle  Park 
North  Carolina  27709 


UNIVERSITY  OF  WISCONSIN— MADISON 
MATHEMATICS  RESEARCH  CENTER 


AN  EXTENSION  OF  THE  KREISS  MATRIX  THEOREM 

Bruce  A.  Wade 
and 

John  C.  Strikwerda 

Technical  Summary  Report  #2951 
August  1986 

ABSTRACT 

A  new  condition  is  shown  to  be  equivalent  to  the  other  conditions  of  the  Kreiss  Matrix 
Theorem  for  power  bounded  families  of  matrices.  This  new  condition  is  important  for  the 
application  of  the  theory  of  pseudo-difference  operators  to  stability  estimates  for  variable 
coefficient  finite  difference  equations.  As  an  example  of  the  usefulness  of  this  new  condition, 
we  use  it  to  prove  stability  of  the  leap  frog  scheme  for  hyperbolic  equations  with  variable 
coefficients. 
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SIGNIFICANCE  AND  EXPLANATION 


The  existence  of  a  symmetrizer  is  a  vital  part  of  the  new  stability  estimate  obtained 
in  !  10]  for  general  variable  coefficient  linear  finite  difference  equations.  It  allows  a  Girding 
inequality  for  pseudo-difference  operators  to  be  utilized  and  is  the  finite  difference  analogue 
of  the  symmetrizer  for  partial  differential  equations.  The  theorem  which  appears  in  this 
paper  links  the  existing  necessary  and  sufficient  conditions  for  the  power  boundedness 
of  a  family  of  matrices  to  the  existence  of  a  matrix  symmetrizer.  An  example  is  given  to 
indicate  how  the  new  condition  allows  stability  estimates  to  be  obtained  for  finite  difference 
schemes  which  do  not  satisfy  the  hypotheses  of  previous  stability  theorems. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies 
with  MRC,  and  not  with  the  authors  of  this  report. 


AN  EXTENSION  OF  THE  KREISS  MATRIX  THEOREM 


Bruce  A.  Wade 
and 

John  C.  Strikwerda 

The  Kreiss  Matrix  Theorem  [2]  gives  necessary  and  sufficient  conditions  for  the  power 
boundedness  of  a  family  of  matrices.  The  power  boundedness  property  is  related  to  the 
stability  of  variable  coefficient  finite  difference  equations  which,  in  turn,  is  related  to  con¬ 
vergence.  In  order  to  obtain  stability  and  convergence  estimates  for  variable  coefficient 
finite  difference  equations  it  is  useful  to  introduce  a  new  condition  in  the  Kreiss  Matrix 
Theorem.  This  condition  involves  the  construction  of  a  matrix  called  a  symmet.rizer. 

The  symmetrizer  is  used  in  jlO]  to  obtain  stability  estimates  for  variable  coefficient 
finite  difference  equations  through  the  use  of  a  Girding  inequality  and  the  theory  of  pseudo¬ 
difference  operators.  A  similar  theory  of  pseudo-difference  operators  for  regular  elliptic 
systems  of  difference  equations  was  developed  by  Bube  and  Strikwerda  in  [l].  Michelson 
;6j  has  also  developed  a  theory  of  pseudo-difference  operators.  The  matrix  constructed  is 
called  a  symmetrizer  because  it  is  the  finite  difference  analogue  of  the  symmetrizer  for  the 
Cauchy  problem  in  the  theory  of  partial  differential  equations.  We  include  an  example  of 
the  application  of  the  new  symmetrizer  condition  to  the  proof  of  the  stability  of  the  leap 
frog  difference  method  for  variable  coefficient  problems. 

Tadmor  9  generalized  the  Lax-Wendroff  5  sufficient  condition  for  stability.  The 
Lax  VNendroff  version  of  the  theorem  involves  the  numerical  radius  of  the  amplification 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041. 


matrix  while  Tadmor's  condition  uses  a  generalized  numerical  radius.  The  generalized 
numerical  radius  condition  was  shown  to  be  necessary  and  sufficient  for  the  power  bound¬ 
edness  of  a  family  of  matrices.  We  show  how  the  new  symmetrizer  condition  is  related  to 
Tadmor’s  condition.  In  the  interest  of  showing  this  relationship,  we  include  more  implica¬ 
tions  in  the  proof  than  are  necessary  to  establish  the  theorem. 

Consider  the  vector  space  Cm  and  let  A  and  fi  be  m  by  m  complex  matrices  with  fl 
Hermitian  and  positive  definite.  Define,  for  x,y  €  Cm, 

(x,y):=x-y,  j'xi2  :=  <*,*). 


and 


||j4;!  sup 


\Ax\ 


x/o  N! 


Define  the  D-numerical  radius  of  A  to  be 


!(n  Ax,x)\ 
rn(A)  :=  sup  — 


j  The  real  part  of  a  matrix  A ,  namely  |  (A'  -  A),  where  »  denotes  the  conjugate  transpose, 

will  be  denoted  by  5JM.  Given  two  Hermitian  matrices  A  and  B,  one  says  that  A  <  B  if, 
for  each  x  €  C"  . 

(A  x,  x)  <  (Bx,x). 

Given  an  infinite  family  of  m  by  m  complex  matrices  denoted  by  J .  the  following  the- 
i  orem  contains  four  known  equivalent  conditions  as  well  as  the  new  symmetrizer  condition. 

l 

\ 

I  The  original  version  was  given  by  Kreiss  in  2  . 
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Theorem  (Kreiss  Matrix  Theorem) 

The  following  are  equivalent: 

A)  There  exists  C a  >  0  such  that  for  all  A  E  J  and  for  each  n  E  Z+,  ||An!|  <  C 
R)  There  exists  Cr  >  0  such  that  for  all  A  E  land  for  all  \z\  >  1, 

W(zl-A)-'i<Cg(\z,-  1)-'. 

H)  There  exists  Cr  >  0  such  that  for  each  A  E  ?  there  is  a  positive  definite,  Hermitian 
matrix  H  which  satisfies  C#1!  <  H  <  CrI  a nd  A’ HA  <  H . 

N)  There  exist  Cr  >  0,  C0  >  0  such  that  for  each  A  E  J  there  is  a  positive  definite, 
Hermitian  matrix  N,  called  a  symmetrizer,  which  satisfies  C^1 1  <  X  <  Cr  and,  for 
each  \z\  >  1. 


X(N{I  -  z~lA))>C0{l-  Izr1)/- 

U)  There  exists  C n  >  0  such  that  for  each  A  E  J  there  exists  a  positive  definite,  Hermitian 
matrix  fl  which  satisfies  Cq1I  <  O  <  CqI  and,  for  all  n  E  Z  +  ,  rn(.4n)  <  1. 

The  original  paper  of  Kreiss  '2'  included  the  conditions  A.  R.  H  and  another  condition 
which  is  commonly  called  S.  We  omit  the  condition  S  for  brevity.  The  condition  H  was 
introduced  by  Tadmor  in  [9|.  The  condition  N,  which  will  be  called  the  symmetrizer 
condition,  is  new.  Although  the  symmetrizer  condition  may  appear  to  be  weaker  that 
Tadmor 's  generalized  numerical  radius  condition,  we  show  directly  that  these  conditions 


are  equivalent,  thus  elucidating  the  relationship  between  them. 


PROOF: 


The  proof  that  condition  A  implies  the  resolvent  condition  R  and  that  R  implies  H  are 
contained  in  ■  2] .  Kreiss  proved  that  R  implies  H  by  first  proving  an  intermediate  condition, 
which  he  called  S.  This  is  the  most  difficult  part  of  the  theorem  and  involves  estimates  on 
the  behavior  of  the  eigenvalues,  especially  for  those  near  the  unit  circle.  Tadmor  proved 
in  [9]  that  fi  implies  A  by  utilizing  the  generalized  Halmos  inequality,  rn(An)  <  (rn(A))n, 
(see  i8,9i),  and  the  inequality  M!  <  2r/(A).  This  last  inequality  can  be  proved  by  writing 
.4  as  the  sum  of  its  Hermitian  and  skew  parts  and  noting  that  the  spectral  radius  of  a 
normal  matrix  equals  its  spectral  norm. 

We  now  prove  that  the  symmetrizer  condition  is  equivalent  to  the  others. 

We  begin  by  prov  ing  that,  condition  H  implies  the  symmetrizer  condition.  N.  Assuming 
H  holds,  let  A  £  7  and  let  H  be  the  corresponding  positive  definite  Hermitian  matrix. 
Expanding  the  relation 

0  <  (1  -  z-lA)'  H  (I  -  z-'A), 

we  have 

0  <  H  -  23?  ( z~'HA )  +  \z  ~2A’HA 
<  H  23 l{z~'HA)  +  ! z\~2H 
-  23?  (H[l  -  z~lA ))  +  (| z\~2  -  1  )H. 

Thus 

I  (1  -  ^  <  |  (l  -  !*!-*)« 

<  D?(//(/  -  z-'A)), 

and  the  symmetrizer  condition  follow's  with  N  taken  as  H . 

The  symmetrizer  condition.  N.  implies  H.  To  see  this,  note  that  if  N  holds  then 

3?  (.\(J  -  z-'A))  >  C0(l  ~  ' z  ~')I 


This  implies  that 


SR  (Ar)  >  SR  (Nz  1A),  for  all  jzi  >  1. 

Given  x  6  Cm,  choose  z  such  that 

SR  ((N Ax,x)z~')  =  j(JV>4 x.x)\. 

Therefore, 

(N A  x,x)\  <  (Nx,x),  for  a]]  x€Cm, 

which  implies  rjv(.4)  <  1.  This,  together  with  the  generalized  Halmos  inequality  [8,9', 
proves  that  condition  fi  holds. 

Conversely,  assume  the  condition  fi  holds.  We  will  show  that  this  implies  the  sym- 
metrizer  condition,  N.  This  part  of  the  proof  is  included  only  to  show  the  relation  between 
the  symmetrizer  condition  and  Tadmor’s  generalized  numerical  radius  condition.  We  have, 
forzeC'". 

5R(fi  (/  -  z~xA)x,x)  =  SR;fix,x)  -  SR(z-1fl.4x,z) 

>  (flx,x)  -  jzj-1|(n.Ax,x' j 

>  (Ox,  x)  -  \z\~x  (Clx.x) 

-  (l- ! srl)ittxtx) 

>  Gn1  (l  -  I*;-1)  ix,x). 

This  implies  the  symmetrizer  condition  with  N  :=  Cl,  and  Cs  :=  C q,  and  Co  :=  Cq1  . 

The  symmetrizer  condition  is  also  closely  related  to  the  resolvent  condition.  If  condi¬ 
tion  N  holds,  then 

C0  (l  -  iz  _1)  <  !|SR  (N(I  -  z~i  A)) 

<  !;ATi :1:/ -  s_,.4  . 

<  C,v i  1  -  z~l A  . 
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where  the  last  inequality  follows  from  observing  J'AT!  =  u/(N)  <  Cs ■  and  the  first  follows 
from  the  Cauchy-Schwarz  inequality.  This  is  easily  seen  to  be  equivalent  to  the  condition 
R.  The  proof  of  the  theorem  is  complete. 

We  now  present  a  simple  example  which  illustrates  the  use  of  the  symmetrizer  condi¬ 
tion  in  proving  the  stability  of  variable  coefficient  finite  difference  equations.  The  methods 
of  Kreiss  [3],  Parlett  ;7j  and  Lax  and  Nirenberg  [4]  do  not  apply  to  the  difference  scheme 
given  here:  however,  the  symmetrizer  condition  can  be  used  to  obtain  the  stability  esti¬ 
mate.  We  make  use  of  the  calculus  of  pseudo-difference  operators  in  the  form  developed 
by  Bube  and  Strikwerda  jl  and  Michelson  [6].  A  complete  presentation  of  the  applica¬ 
tion  of  the  symmetrizer  condition  to  stability  and  convergence  estimates  of  general  finite 
difference  equations  along  with  applications  will  be  presented  in  [10  . 

Consider  the  variable  coefficient  initial  value  problem 

ut  ~  a(x,t)ux ,  t  >  0,  j  £  R 

where  uo(-)  £  L2(R)  is  a  scalar  valued  function  and  a(x,t)  is  constant  outside  of  a  compact 
set.  We  approximate  the  partial  differential  equation  by  the  leap  frog  difference  scheme 


rB+I  -  vn'  1 

?>i  m 

“Wn  fcl  rm-l)t 

n  >  1,  m  6  Z 

Vm  : 

=  u0(mh) 

(1) 

v'm: 

ft  1  <  /  o  0 

=  Vm  +  2^am.n  (*W  1  ~  Vm- 

i)  * 

where  a m.n  a(mh.  nk)  and  A  k/h  is  constant.  We  assume  that  sup,.  #  A!a(x,t)l  <  1. 

Define  the  discrete  function 

f  (l/2fc)v£,.  if  n  =  -1; 

A”  :=  -j  (1/2*)  (»“  -  -  r^.,)).  if  n  =  0; 

l  0.  else, 
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then  it  is  easy  to  check  that  the  difference  equation  (1)  is  equivalent  to 


,n+  1 


-  Vr 


—  ^ am,n  (t* 


tn-t-  1 


C-,) 


n,m  e  Z, 


(2) 


where  we  now  consider  the  difference  equation  as  being  for  all  time  levels  and  take  vn  to 
be  0  if  n  is  negative.  Note  that  essentially  acts  as  a  “delta  function”and  its  nonzero 
elements  are  of  order  1/A:. 

We  now  reduce  this  difference  scheme  to  a  two  level  scheme.  If  we  define 


and 


equation  (2)  may  be  written  as 


W »  -  G(m, «)»'*- 1  +  kF*.  n,m  €  Z 


(3) 


where 

G(m,n)~  ^.»(T-T-')  lj 

and  T  is  the  translation  operator  in  x. 

The  symbol  of  the  difference  operator,  G ,  is  defined  to  be 


G(m,n,Z) 


(  2Atam  „  sin  $h.  1  \ 

V  1  o)' 


We  note  that  for  each  value  of  (m,n)  the  eigenvalues  of  G(m.n,  £)  are  on  the  unit  circle,  in 
particular  the  scheme  is  nondissipative.  Therefore,  the  method  of  Kreiss  3  and  Parlett  1 7 ; 
for  proving  stability  does  not  apply.  The  sharp  Girding  inequality  of  Lax  and  Nirenberg 


>"C> 


[41  also  does  not  apply  since  the  norm  of  the  symbol  is  not  bounded  by  1.  We  will  utilize 
the  symmetrizer  condition  in  order  to  prove  that  this  difference  scheme  is  stable. 

The  theory  of  pseudo-difference  operators  is  employed  in  the  context  of  the  discrete 
Sobolev  space 

HVth  :=  {v  :  hZ  x  kZ  ->  R  :  6  L2(Dr?i/l)}, 

where  D,,k  :=  {(£.r)  :  if!  <  |r|  <  f },  s  =  r;  -t  ir,  with  rj  >  ij0  >  0,  and  t>(f,s) 

denotes  the  discrete  Fourier  transform  in  space  and  Laplace  transform  in  time  of  defined 
by 


v  (£,«)  = 


hk 

2ir 


V''  .-ank-tfmA.n 


n.m—  -  oe 


The  corresponding  norm  is 


=  hk  \e 


m,n  =  -oo 


See  Michelson  ;6i  or  Wade  10;  for  details  on  the  pseudo-difference  operator  theory  as  well 
as  the  discrete  transform  and  discrete  Sobolev  space. 

Consider  the  family  of  matrices 


5 

V 


7  {e  rt,,kG(m,n,  £)  :  m,n  c  Z,  :£i  <  -}, 

where  r?0  is  any  positive  number.  Since  the  eigenvalues  of  every  element  of  7  are  distinct, 
inside  the  unit  circle  and  uniformly  separated,  it  is  clear  that  the  resolvent  equation  holds. 
Therefore,  there  is  a  constant  C\,  >  0  such  that  for  each  element  of  7  there  is  a  symmetrizer. 
S(m.  n,  £).  which  satisfies 


3?  ( N{m.  n.  £ )(/  -  e  pkG(m.n.£)^J  -e  vk )  1 .  rj  >  rj0  >  0. 


(4) 
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We  now  apply  a  Girding  inequality  for  pseudo-difference  operators,  (see  ;lj,  j6i,  or 
[  10]) .  The  Girding  inequality  essentially  says  that  the  condition  (4),  which  holds  for  the 
symbol  of  the  difference  operator,  is  also  true  for  the  variable  coefficient  operator  modulo 
lower  order  terms.  We  have 


1  -  e~r,k 


)  WWU  <  JV(J  -  G(n,m)S~l)W)^h  -r  c  \W |!*.fc  (5) 


where  5  is  the  shift  operator  in  t. 

By  choosing  r?o  larger,  and  the  constant  Co  smaller,  and  using  that  W  satisfies  the 
equation  (3),  this  implies 


Co  * — l - WWU  <  *WKF)n,h  <  c  w°\l  V  >  Vo- 


The  last  inequality  follows  from  the  fact  that  the  grid  function  F  is  supported  on  only  two 
time  levels.  In  terms  of  the  original  variable  r"  this  is 


1  -  €~,lk  '' 


)  MliU  £  ci>  T. 


r?  >  J?o  >  0, 


for  some  constant  C. 


The  estimate  (6)  is  equivalent  to  the  standard  notion  of  stability.  By  ‘standard  notion 
of  stability'  we  mean  the  definition  which  requires  the  space  norm  of  the  solution  at  any 
time  level  n  to  be  bounded  by  a  constant  depending  on  nk  times  the  space  norm  of  the 


initial  data  of  the  difference  scheme.  The  reader  is  referred  to  the  paper  of  Bube  and 


Strikwerda  1  for  a  similar  estimate  which  uses  a  G&rding's  inequality  for  elliptic  systems. 

Finally  we  comment  on  the  reason  that  the  symmetrizer  condition  allows  estimates 
for  more  general  difference  schemes.  The  previous  methods  for  proving  stability  involved 


-.v.% 


Wl 


ST*.' 


1 


$2 


w 


•’.VI 


*.*!‘.*l 


the  conditions  H  or  0.  These  conditions  involve  nonlinear  operations  on  the  elements  of 


the  family  of  matrices,  while  condition  N  is  essentially  a  linear  operation  on  the  elements 


of  the  family.  The  condition  H  contains  the  term  A’  HA  and  the  condition  fi  contains 


the  absolute  value  of  (H*,-),  both  of  which  are  nonlinear  operations  on  the  matrices  in  the 


family.  The  symmetrizer  condition  utilizes  only  the  multiplication  of  elements  of  the  family 


by  a  matrix,  in  addition  to  the  operation  of  taking  the  real  part  of  a  matrix.  The  use  of 


only  linear  relations  allows  for  greater  flexibility  in  the  application  of  Girding  inequalities. 


Estimates  for  more  general  schemes  as  well  as  other  results  are  contained  in  j  10]. 


i  » 
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